co 

o 

(N 

S-H 

in 



s 



On the Griinwald-Letnikov fractional-order derivative 
with fixed memory length 

Mohammed-Salah Abdelouahab*, Nasr-Eddine Hamri 

Department of Mathematics and Computer sciences 
Mila University Centre, 43000, Algeria 



Abstract 



Contrary to integer order derivative, the fractional-order derivative of a non- 
constant periodic function is not a periodic function with the same period, 



on 
Q 

as a consequence of this property the time-invariant fractional order systems 



does not have any non-constant periodic solution unless the lower terminal 



of the derivative is ±00, which is not practical. This property limits the 

,_i applicability areas of fractional derivatives and makes it unfavorable, for a 

wide range of periodic real phenomena. Therefore enlarging the applicability 

of fractional system to such real area is an important research topic. In 
CO 

this paper we give a solution for the above problem by imposing a simple 
# > 

modification on the Grunwald-Letnikov definition of fractional derivative, 

s_i 

this modification consists of fixing the memory length and varying the lower 

terminal of the derivative. It is shown that the new proposed definition of 

fractional derivative preserves the periodicity. 
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1. Introduction 

In 1695, Leibniz laid the foundations of fractional-order derivative which 
means the extension of integer-order derivative concept pQ, but the first 
systematic studies seem to have been made at the beginning and middle 
of the nineteenth century by Liouville, Riemann, and Holmgren. Liouville 
has expanded functions in series of exponentials and defined the nth-order 
derivative of such a series by operating term-by-term as though n were a 
positive integer. Riemann proposed a different definition that involved a 
definite integral and was applicable to power series with non-integer expo- 
nents. It was Griinwald and Krug who first unified the results of Liouville 
and Riemann. Griinwald, by returning to the original sources and adopting 
as starting point the definition of a derivative as the limit of a difference quo- 
tient and arriving at definite-integral formulas for the nth-order derivative. 
Krug, working through Cauchy's integral formula for ordinary derivatives, 
showed that Riemann's definite integral had to be interpreted as having a 
finite lower limit while Liouville's definition corresponded to a lower limit 
— oo [2H6]. It turns out that the Riemann-Liouville derivatives have certain 
disadvantages when trying to model real-world phenomena with fractional 
differential equations. In the second half of the twentieth century Caputo 
proposed a modified concept of a fractional derivative [7], by inverting the 
composite order of integer-order derivative operator and integral operator 



arising in the Riemann-Liouville definition, the resulting definition is better 
suited to such tasks. It has been found that many systems in interdisciplinary 
fields can be described by the fractional differential equations, such as vis- 
coelastic systems, dielectric polarization, electrode-electrolyte polarization, 
electromagnetic waves, and quantum evolution of complex systems [8HT2]. 
Due to the growing interest of fractional-order derivatives to be applied in 
different areas, it seems analysis of this type of derivative is of great im- 
portance. The existence of periodic solutions is often a desired property in 
dynamical systems, constituting one of the most important research direc- 
tions in the theory of dynamical systems, with applications ranging from 
celestial mechanics to biology and finance. In chaos theory, the idea of chaos 
control is based on the fact that chaotic attractors have a skeleton made 
of an infinite number of unstable periodic orbits which are subjected for 
stabilization [T3]. Thus for chaos control in fractional order systems, it is 
important to show that periodic solutions exist in such systems. Tavazoei 
[H] have proved that the fractional-order derivatives of a periodic function 
with a specific period cannot be a periodic function with the same period, as 
a consequence, the periodic solution cannot be detected in fractional-order 
systems, under any circumstances [T5T - I2C)] . This property limits the applica- 
bility areas of fractional-order systems and makes it unfavorable, for a wide 
range of periodic real phenomena. Therefore enlarging the applicability of 
fractional-order systems to such real area is an important research topic. The 
most important contribution of this paper is to this end. 



2. Griinwald-Letnikov fractional derivative 

Lets a > 0, the Griinwald-Letnikov ath order fractional derivative of 
function /(£) with respect to t and the terminal value a is given by [2] 
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where T is the gamma function defined by the Euler limit expression 



T(x) = lim 
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n-s-oo x(x + l)...(x + n) 

where x > 0. 

or the so-called Euler integral definition: 
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The following theorem gives negative answer for question of preservation 
of the periodicity 



Theorem 1. Suppose that f(t) is (m — l)-times continuously differentiable 
and f( m > is bounded. If f(t) is a non-constant periodic function with period 
T,then the functions G ^D^f(t), where < a ^ N and m is the first integer 
greater than a, cannot be periodic functions with period T . 

Example 2. The derivative of the sine function is given by 

G L D?sin(t)=t 1 - a E 2t2 _ a (-t 2 ), 

where E a ^ is the two-parameter Mittag-Laffter function defined by /||/ 

°° t k 

Figure HI illustrate numerical approximation of ^Df^sin^t) which is not 
periodic but it converges to the periodic function sin(t + a|). 

As a consequence of the above theorem, the periodic solution cannot be 
detected in fractional-order systems, under any circumstances 



Corollary 3. A differential equation of fractional-order in the form 

G a L D?x(t) = f(x(t)), 
where < a ^ N, cannot have any non-constant smooth periodic solution. 



This property makes fractional systems not applicable for a wide range 
of real periodic phenomena. 

3. Griinwald-Letnikov fractional derivative with fixed memory length 

In this section we introduce a new definition of fractional derivative based 
on the Griinwald-Litnikov definition and we call it the Griinwald-Letnikov 
fractional derivative with fixed memory length. 

Definition 4. Let a > 0, L > 0, m an integer such that m — 1 < a < m and 
f an integrable function in the interval [a — L,b]. We define the Griinwald- 
Letnikov fractional derivative with fixed memory length L of f , Dff by 
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D?f(t) = lim -L Y"(-l) fc , IT 5 (a + , 1 \, f(t - kh), te[a,b]. (5) 
tJKJ h ^oh aZ ^ K ' k\T(a - k + 1) v ' L J v; 



Remark 5. A remarkable difference between classical definitions of fractional- 
order derivative and the new definition is that in the classical definitions the 
lower terminal a is fixed and the memory length t — a is a function of the 
variable t but in the new definitions the lower terminal a = t — L is a function 
of the variable t and the memory length L is fixed. 

The definition of the Griinwald-Letnikov fractional derivative with fixed 
memory length L of /, given in definition Q is expressed in term of a limit 
which is not suitable for analytic study. The following proposition gives an 
evaluation of this limit under some restrictions of the function /. 
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Proposition 6. Under the assumptions of definition H\) if the function f is 
(m + l)-differentiable with /( m+1 ) e L\[a — L, b] we have 
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where we denote 



A/(t - fc/i) = /(t - JWi) - f(t -{k + l)h). 



Applying the property ( 10 ) of binomial coefficients repeated m times, starting 



from (13) we obtain: 
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Let us evaluate the limit of the fcth term in the sum ( 14 ) 
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Taking into account (17) and (18) and applying theorem (2.1) in [2] we obtain 
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Proposition 7. (Fractional derivative of a power function) 

Let a > 0, L > 0, m an integer such that m — 1 < a < m and f(t) = t n , 

then we have 
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Proof. If n < m, then /( m+1 ) = 0, substituting in won yields the relation (20). 



If n > m, then from (pi) and using (n — m — 1) successive integration by part 



we obtain the relation (20). 

□ 

Corollary 8. (Fractional derivative of a constant function) 

If f is a constant function (i.e. f(t) = C for allt G [a — L,b}), then iDff 

is a constant function for all t G [a, b]). Furthermore we have 

"w-itw^r (21) 

Proof. Suppose that / is a constant function ( f(t) = C for all t G [a — L, b]) 
, then f k \t-L) = /^ m+1) (r) = for k = 1,2, ..., m substituting in (§) yields 
the relation (1211) □ 



Proposition 9. (Fractional derivative of the exponential function) 

Let a > 0, L > and m an integer such that m — 1 < a < m, then we have 



lAV) = %^. (22) 
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3.1. An interpolation property 

The fractional-order derivative mean the extension of the integer-order 

derivative to real order, in other word the fractional-order operator D? is 

d n 
an interpolation of the integer-order operator — , so it is naturally to ask 

dt n 

the question: is the proposed Griinwald-Letnikov fractional-order derivative 
with fixed memory length an extension of the integer-order derivative? The 
following proposition gives the answer to this question. 

Proposition 10. Let a > 0, L > 0, m an integer such that m — 1 < a < m 
and f an m-differentiable function on [a — L, b]. Then for all t G [a, b] we 
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Thus the fractional-order operator L Df is an interpolation of the integer- 
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d n 
order operator - — . D 

F dt n 

Example 11. Let a > 0, L > 0, m an integer such that m — 1 < a < m 

and suppose that fit) = C is a constant function then from corollary M) we 

have 

lim M lD?f(t) = lim j^ r = = /<™)(t), (25) 



and 



lim "lD ? f( t ) = Inn ] -^— - = C = /(«). (26) 



5.^. Derivative of periodic function 

The main result of this paper as previously announced is the introduction 
of new modified definition of fractional derivative with finite fixed memory 
length which preserve the periodicity, this result is summarized in the fol- 
lowing theorem 

Theorem 12. Let a > 0, L > 0, m an integer such that m — 1 < a < m 
and f a function defined on [a — L, b] such that tP'tf exists. If f is a 
periodic function with period T (i.e. f(t + T) = f(t) for all t e [a — L, b]), 
then jfltf ^ s a P e ™dic function with the same period T. 

Proof. Suppose that / is a periodic function with period T, then we have 
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Hence i^tf * s a periodic function with period T. 

□ 



Example 13. (Fractional derivative of sinus and cosines functions) 

oo £2p+l 

1) We have sin(t) = ^ ( — l) p — — -, then 

p =o (2p + 1)! 



oo MG 

MG ^„ . , , V"^ 



L * U ^ ^ (2p + l)! 
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using the relation (20) we obtain 



oo 2p+l rk-a(f_ r\2p+l-k 

L D?sin(t) = £ E (~l) p r9 ^ l ' /, TTT 
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r(4-a) 

oo J_^V- a oo ^2p+l-Q 

sm(i - L) V (-l) p — r + cosit - L) Y\ (-l) p — r 

= L- a sin(t - L)E 2 ^ a (-L 2 ) + L l - a cos{t - L)E 2 ^ a (-L 2 ) 
A.sin(t - L) + B.cos(t - L), 



where A = L~ a E 2 i_ a (-L 2 ) and B = L y - a E 2 2 J-L 



MO 



Clearly L Dfsin(t) is a periodic function with period 2tt. 

Figure \!2\ shows the modified fractional- order derivative L D^sin{t) of 
sine function for L = 30 and some value of a which confirms the analytic 
result obtained. 

Remark 14. When the memory length L increases its influence on the 
derivative decreases, as shown in Figure [3| thus in order to obtain an in- 
dependent fractional derivative from L it is recommended to chose L greater 
as possible. 
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2) We have cos(t) = £(-1) p t^t7, then 
p =o (2p)' 



P =o (2p)! 

oo 2p L k- a U_ L Y2 V ~k 

- E E(-i) p ^ 



P =ofc=o (2p-k)\T(k-a + l) 

j —a rl-a j 2— a 
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r(l-a) v ' r(2-a) v ' r(3-a) 

+ f[4^a) Sin{t - L) + - 

oo j2p-a oo j2p+l-a 

= cos(t -L)f] {-l) p — r - sin(t -L)T (-l) p — r 

= L~ a cos{t - L)E 2>1 _ a (-L 2 ) + L l - a sm{t - L)£ 2 , 2 _ Q (-L 2 ) 
-- A.cos{t - L) - B.sin(t - L). 



Clearly L D"cos(t) is a periodic function with period 2-ir. 






Figure 1: Fractional derivative of sine function Dfsinif) for a = 0.5 



y= Z **? sin '») 




Figure 2: Fractional derivative of sine function L Df sin{t) for L = 30 and some value 
of a. 



y= ™'D' t J " sin(t> 




Figure 3: Influence of the memory length L on the fractional derivative L D^ 5 sin(t) 
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4. Conclusion and discussion 

Although the idea of fixed memory length is inspired from the short- 
memory principle introduced by I.Podlubny (1999) for numerical needs, there 
is a significantly difference between this two ideas near the starting point. 
Namely for the calculation of the fractional derivative of a certain function /, 
on the interval [a, b], using the first idea, the memory length L is fixed for all 
t G [a, b], as previously mentioned, but using the second idea, if t G [a,a + L] 
then the lower terminal is fixed at a, and the memory length L is a function 
of t ( L = t — a ) and this do not allow the preservation of periodicity. The 
modified definition of Griinwald-Letnikov fractional-order derivative reported 
in this paper posses two useful properties, the first is the preservation of 
periodicity as it is demonstrated and the second is the short memory which 
reduces considerably the cost of numerical computations. In the forth coming 
papers we extend this modification to the other definitions of fractional-order 
derivative (Riemann-Liouville, Caputo...). Now, this related question may 
be raised. Can autonomous fractional differential systems with fixed memory 
length generate exactly periodic signals? Finding the answer of this question 
and investigate other property of the proposed fractional derivative, can be 
an interesting topic for future research works. 
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